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ABSTRACT 

A  Random  Network  model  Is  Introduced  to  simulate  the  porous  metal- 
electrolyte  interface.  This  is  the  first  time  that  an  attempt  has  been 
made  at  defining  a  model  for  this  system  that  maintains  the  random 
nature  of  the  rough  topology.  A  previously  utilized  model  for  this 
system  -  the  single  pore  model  is  explored  and  extended.  The  models 
are  compared,  and  the  results  of  the  Random  Network  model  are  found 
to  be  in  qualitative  agreement  with  the  single  pore  models. 
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INTRODUCTION 


The  porous  metal  electrode  is  a  subject  of  intense  theoretical  and 
applied  interest.  Although  widely  utilized  in  electrochemical  systems, 
especially  in  battery  technology,  very  little  is  known  in  detail  about  the 
porous  interface  from  first  principles,  aside  from  a  general  knowledge 
of  how  the  porous  electrode's  behavior  deviates  from  flat-electrode 
behavior.  Specifically,  it  would  be  useful  to  know  how  the  complicated 
topology  of  the  porous  metal  electrode  contributes  to  the  electrical 
properties  of  the  interface.  Previous  attempts  at  modeling  the  porous 
electrode  are  presented  in  the  review  article  by  DeLevie  [1].  The 
macroscopic  model  [2]  treats  the  elctrode -electrolyte  system  as  a 
superposition  of  two  continue,  one  of  the  electrode  matrix  and  one  of 
the  solution  matrix  that  fills  all  the  unoccupied  space.  The  solution  is 
assigned  an  effective  conductivity,  and  the  interface  is  considered  to 
have  an  effective  capacitance  per  unit  area.  The  single  pore  model  is 
based  on  a  one-dimensional  representation  of  a  pore  as  a  transmission 
line  (see  figure  S)  which  represents  a  pore  of  "average  pore  length". 
These  models,  however,  do  not  account  for  the  random  distribution  of 
pores  in  a  true  porous  solid,  nor  for  a  complicated  surface  topology. 
We  have  chosen  to  introduce  a  three  dimensional  model  of  the  rough 
Interface  based  on  the  random-network  lattice  that  has  proven  to  be  so 
useful  in  representing  conduction  in  Inhomogeneous  materials  and 
related  percolation  problems  [3,4,5],  A.c.  impedance  measurements  are 
widely  used  experimentally  for  the  in -situ  characterization  of 
electrochemical  systems  under  equilibrium  conditions.  In  conducting 
these  measurements  on  porous  electrodes,  the  complex  topology  of  the 
electrode  in  various  states  of  charge  yields  ambiguous  data,  not  subject 
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to  a  unique  interpretation.  The  computer  simulation  allows  us  to 
examine  the  impedance  of  the  system  as  a  function  of  given  charge 
accumulation  modes  at  the  electrode -electrolyte  interface.  The  model  is 
not  intended  to  portray  the  d.c.  operation  of  the  electrode,  but  rather 
to  provide  a  means  for  understanding  the  interface.  Using  an  IBM/370 
computer,  we  construct  our  model  and  calculate  its  complex  impedance, 
Z(w),  in  the  frequency  range  of  one  Hertz  to  one  Megahertz.  Based  on 
this  impedance  spectrum,  we  model  the  system  in  terms  of  passive 
elements  [6] . 

In  this  paper,  we  will  present  our  model  in  detail,  discuss 
preliminary  results,  and  compare  these  results  with  the  "DeLevie 
model”.  We  will  also  discuss  further  work  being  done  to  utilize  this 
model  for  a  more  detailed  characterization  of  the  interface. 
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I.  THE  "RANDOM"  NETWORK  MODEL 
We  consider  a  simple  system,  in  which  a  1  cm3  porous  metal 
electrode  is  immersed  in  electrolyte,  along  with  a  counter* electrode. 
Only  the  bottom  (1  cm1)  surface  of  the  electrode  is  exposed  to  the 
electrolyte.  The  electrode  is  allowed  to  discharge,  allowing  for  a  build¬ 
up  of  a  semiconductor  on  the  metal  surfaces  exposed  to  the  elecrolyte. 
To  explore  the  electrical  behavior  of  this  system,  we  note  that  we  may 
define  a  local  conductivity,  o(f),  the  bulk  conductivity  at  point  r  in 
our  system.  Rather  than  attempt  to  solve  for  the  current  distribution 
using  the  electrostatic  equations  for  the  continuous  case,  we  transform 
the  problem  into  one  requiring  the  solution  of  a  discrete  set  of 
Kirchoff's  Law  equations  (a  finite  difference  approximation  to  the 
continuum  problem)  [3] .  We  subdivide  our  bulk  system  into  a  regular 
cubic  mesh  of  points  {i^}  (1=1,  S)  and  assign  to  each  branch  of  the  mesh 
a  conductivity  that  represents  the  conductivity  of  the  surrounding  bulk 
in  that  region  of  space.  Let  the  conductance  of  the  circuit  branch  that 
connects  node  i  and  j  be  denoted  as  gy,  the  voltage  at  node  i  be  Vj, 
and  the  net  current  into  node  i  be  Ij.  We  then  have  (S-l)  equations  of 
the  form: 


li*u  <vrV  ■  xi  « 

or  in  matrix  notation: 

G  •  V  *  I 

where  G  is  the  conductance  matrix  (Gy  =  Ij  gy  and  Gy  =  -gy),  V  is 
a  vector  containing  the  voltages  on  each  node  with  respect  to  ground, 
and  1  is  a  vector  containing  the  net  current  Inputs  to  each  node.  We 
may  solve  these  simultaneous  equations  for  the  V's  on  all  the  nodes. 


and  for  the  total  impedance  of  the  ayatem.  A  detailed  presentation  of 
the  algorithm  that  we  have  employed  follows. 

We  define  a  three  dimensional  cubic  lattice,  A(i,j,k),  where  i, J  and 
k  range  from  1  to  N,  and  N3=S,  the  total  number  of  nodes.  N  is 
chosen  so  that  it  is  sufficiently  large  to  approximate  an  infinite  system. 
Construction  of  the  electrode  Is  simulated  using  a  random  number 
generator  to  place  "metal"  at  various  points  in  the  lattice,  until  the 
desired  porosity  is  reached.  A  check  is  made  of  the  resulting  electrode 
to  insure  that  it  is  continuous  and  that  no  pieces  of  metal  are 
"hanging"  in  mid-air.  Any  disjointed  pieces  of  metal  are  removed,  and 
are  randomly  replaced  on  the  remaining  available  lattice  points.  This 
process  is  repeated  until  a  continuous  electrode  of  the  desired  porosity 
is  obtained.  Lattice  plane  k=l  is  defined  as  electrolyte,  to  represent 
the  solution;  the  electrode  itself  begins  at  lattice  plane  k=2.  The 
electrolyte  is  allowed  to  "seep  into"  the  pores  in  the  electrode  until  all 
accessible  pores  are  filled,  by  searching  for  continuous  electrolyte 
paths  from  the  k=l  plane  throughout  the  electrode.  Any  point  on  the 
lattice  which  has  been  left  empty  (i.e.  no  metal  was  placed  there,  and 
electrolyte  was  precluded  from  penetrating  to  that  site)  is  considered  to 
be  "air". 

Each  lattice  point  represents  a  node  in  a  three  dimensional  circuit 
network  (See  figure  1),  and  each  pair  of  nodes  defines  a  circuit  branch 
in  the  network.  The  impedance  of  each  branch  is  determined  by  the 
characteristics  of  the  two  surrounding  nodes.  Thus,  for  example,  if 
A(l,l,l)  is  electrolyte  and  A(l,2,l)  is  also  electrolyte,  the  circuit 
branch  defined  by  those  nodes  will  consist  of  two  series  resistors  of 
resistance  Rfl  (the  resistance  of  a  microscopic  section  of  the 
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electrolyte).  If  A(l,l,2)  is  metal  and  A(l,2,2)  is  metal,  then  the 

circuit  branch  defined  by  those  two  nodes  will  be  two  series  resistors 

of  value  (the  resistance  of  a  microscopic  section  of  the  metal) . 

Finally,  the  branch  defined  by  A(l,l,l)  -  electrolyte,  and  A(l,l,2)  - 

metal,  is  represented  by  an  R  and  an  R  resistor  in  series,  with  the 

m  e 

addition  of  a  parallel  R-C  element  in  series  with  them  to  represent  the 
semiconductor-electrolyte  interface  that  results  (see  Figure  2).  Any 
circuit  branches  leading  into  nodes  that  are  defined  as  "air”  are  taken 
to  be  of  infinite  resistance. 

The  values  that  we  have  chosen  to  use  as  our  unit  impedances  for 
the  individual  components  of  metal,  electrolyte,  and  semiconductor  are 
presented  in  Figure  2.  These  values  were  obtained  by  taking  the 
values  of  bulk  impedance  for  a  1  cm3  sized  sample,  (we  chose  numbers 
characteristic  of  Zn  and  ZnO)  and  scaling  them  down  to  the  magnitude 
of  a  unit  pore  size-  approximately  IOji  in  diameter.  We  define  a 
characteristic  time  constant  of  the  interface,  t  =  RpC.  These 
components  are  assembled  into  a  system,  as  described  above  (see  figure 
3)  and  the  total  impedance  of  the  system  is  calculated.  The  total 
impedance  is  then  rescaled  by  an  appropriate  factor  that  normalizes  the 
whole  system  to  a  1  cm*  size.  Thus,  an  attempt  is  made  to  obtain 
impedance  data  with  numbers  representative  of  that  which  might  be 
encountered  in  a  real  system. 

Once  the  electrical  network  has  been  defined,  the  computer  sets  up 
the  Kirchoff's  law  equations.  Due  to  the  fact  that  our  model  includes 
reactive  circuit  elements,  the  quantities  G,  V  ,  and  I  In  equations  (I) 
are  complex,  and  may  be  represented  as  2S-1  equations  of  the  form: 


l  UWjUWBU 


(Gr  -Gi  I  I  Vr  I  I  Ir  | 

Gi  Gr  J  l  Vi  J  =  Ini  (2) 

where  the  subscripts  r  and  1  stand  for  the  real  and  Imaginary 
components  respectively. 

A  known  current  is  sent  uniformly  into  all  the  bottom  nodes  (so  as 
to  eliminate  the  "edge  effects"  that  would  result  if  the  current  were 
sent  only  into  one  node),  and  the  impedance  of  the  system  is  calculated 
between  one  node  on  the  electrolyte  plane  (1,1,1)  and  one  node  of  the 
top  electrode  plane  (N,N,N)  by  solving  for  the  V's  on  all  the  nodes. 

For  a  system  of  N=15  (a  lattice  of  size  15x15x15  with  3375  nodes) 
there  are  up  to  6748  equations  to  be  solved.  This  is  a  formidable 
computer  problem.  Indeed,  the  task  would  be  prohibitive  without 
taking  advantage  of  the  fact  that  the  conductance  matrix,  G,  is  a 
sparse-symmetric  coefficient  matrix. 

A  number  of  techniques  exist  for  the  solution  of  this  class  of 
problems.  A  widely  used  method  for  solving  the  simultaneous  equations 
in  the  Random  Resistor  lattice  is  the  Gauss  Seidel  iteration  procedure 
with  over-relaxation  [5] .  This  is  a  very  economical  method,  in  both 
execution  time  and  storage  space.  A  serious  drawback  of  this 
technique  is  the  relatively  strict  requirements  that  it  places  on  the 
coefficient  matrix  in  order  for  convergence  to  be  achieved,  namely  that 
the  matrix  be  either  diagonally  dominant,  or  at  least  positive  definite 
[7].  For  Kirchoff's  law  problems  with  purely  resistive  components, 
diagonal  dominance  is  assured,  since  the  diagonal  elements  of  the 
conductance  matrix  are  simply  the  sum  of  the  off  diagonal  elements .  In 
our  problem,  with  complex  Impedances,  we  find  that  at  frequencies  for 
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which  the  real  conductance  approaches  the  same  order  of  magnitude  as 
the  imaginary  conductance,  the  Gauss -Seidel  procedure  does  not 
converge.  This  is  due  to  the  form  of  the  equation  in  (2)  where  there 
are  many  more  off  diagonal  elements  due  to  the  Gj  values.  When  these 
Imaginary  components  are  large,  the  matrix  is  no  longer  diagonally 
dominant,  nor  positive  definite. 

For  most  of  our  computations,  we  have  resorted  to  the  use  of  a 
Gaussian  elimination  routine  which  takes  advantage  of  both  the 
sparseness  and  the  symmetry  of  the  G  coefficient  matrix  [6] .  Only 
the  nonzero  elements  In  the  upper  half  triangle  of  the  matrix  are 
stored,  and  an  efficient  pivoting  strategy  is  chosen  to  minimize  non¬ 
zero  matrix  fill  during  the  pivoting  and  to  minimize  the  number  of 
multiplications  required  in  the  solution.  In  practice,  the  technique  is 
approximately  an  order  of  magnitude  more  expensive  to  use  than  Gauss  - 
Seidel  in  both  speed  and  storage  requirements.  The  advantage  of  using 
the  Gaussian  elimination  technique  is  that  a  solution  is  guaranteed  for 
almost  any  problem.  The  solutions  obtained  from  the  Gauss- Seidel 
procedure  (in  the  range  in  which  convergence  is  reached)  are  in 
complete  agreement  with  the  results  obtained  using  the  modified 
Gaussian  elimination  technique. 

Figure  3  presents  a  system  constructed  with  a  porosity  of  0.30,  and 
the  resulting  impedance  of  the  system  as  a  function  of  frequency  is 
plotted  in  figure  4.  We  have  identified  three  basic  regions  of  interest 
in  the  frequency  range  spanned  by  our  results.  Region  I  is  the  low 
frequency  range,  in  which  our  system  seems  to  behave  as  a  parallel  R- 
C  combination  (See  figure  4).  Region  II  is  in  the  mid-frequency  range, 
with  an  impedance  spectrum  that  cannot  be  represented  using  passive 
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elements,  and  Region  III  is  representative  of  a  series  R-C  combination. 
The  variation  of  the  impedance  data  as  a  function  of  system  size  is 
presented  in  Figure  5  from  size  N=5  to  N=15.  For  the  results  presented 
here,  N=15  was  used  to  ensure  proper  statistical  representation  of  an 
infinite  system. 

II.  THE  SINGLE- PORE  MODEL 

It  is  interesting  to  compare,  at  least  qualitatively,  the  results  from 
our  three  dimensional  model  with  the  results  obtainable  from  the  single 
pore  model  previously  mentioned  [1] . 

De  Levie,  following  on  the  heels  of  other  researchers  (Daniel' -Beck, 
[9]  et.  al.),  suggested  that  each  pore  in  a  porous  electrode  may  be 
thought  of  as  having  a  uniformly  distributed  electrode  and  electrolyte 
resistance  throughout  its  length.  As  such,  the  single  pore  may  be 
modeled  as  in  figure  6a,  where  R  is  the  resistance  per  unit  length  of 
the  electrolyte  solution  inside  the  pore,  and  Z  is  the  impedance  of  the 
electrode-electrolyte  interface.  The  current-voltage  relationships  of 
this  circuit  may  be  expressed  as  a  differential  equation  and  a  solution 
may  be  obtained  for  the  total  impedance  of  the  circuit  as  [1] : 

Zo  =  (RZ)*  coth  (pL) 

where  R  and  Z  are  defined  above,  p=(R/Z)^,  and  L  is  the  length  of  the 
pore.  This  is  the  main,  and  currently  utilized  [10],  result  from  this 
model.  The  main  characteristics  of  the  calculation  are:  (1)  Although 
flat  electrode  impedances  vary  as  a  fuction  of  Z,  porous  surface 
impedances  will  be  dependant  on  (RZ)^ ,  and  (2)  the  contribution  of 
the  electrode  surface  deep  inside  the  pore  is  negligible,  1/p  becoming 
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the  characteristic  "penetration  depth*.  If  pL  »  1  the  pore  behaves 
like  a  semi-infinite  one. 

In  order  to  utilize  these  results  for  comparison  with  the  Random 
Network  Model,  we  must  introduce  a  value  for  Z,  and  we  must  also 
combine  the  single  pore  transmission-line  impedances  to  form  a  complete 
rough  electrode.  The  former  is  accomplished  by  defining  the  electrode  - 
electrolyte  interface  as  a  standard  parallel  R-C  model,  as  depicted  in 
figure  6b.  We  have  used  the  same  numerical  values  for  the  individual 
components  as  those  in  the  Random  Network  Model  (here  they  are  shown 
already  normalized  to  the  lOp  single  pore  size).  As  a  result, 


2  =  rs  +  :  2i 

1  +  OJ  T 
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A  log-log  plot  of  the  relationship  between  Z  and  f  (for  a  single 
pore)  Is  presented  in  figure  7.  It  is  Interesting  to  note  the  similarity 
in  shape  between  the  impedance  spectrum  of  the  single -pore  model 
(even  though  we  are  only  considering  a  single  pore)  and  that  of  our 
Random  Network  Model  (figure  4) .  The  advantage  of  considering  a 
single  pore  at  the  present  stage  is  that  it  is  the  only  system  for  which 
we  have  an  analytic  result  (equation  5),  which  may  be  used  to  derive 
expressions  for  the  impedance  in  different  frequency  ranges. 

We  define  three  frequency  regimes: 

(a)  The  low  frequency  range  (wt  «  1): 

Re(Zo)  »  K  (7a) 

Im(Zo)  *  Kux 


where  K  (independent  of  frequency)  is  given  as: 


e2x  _  e-2x 
e2x  +  e~2x  +  2 


and  x  is: 

x  «  L  (R/Rp)* 


The  notable  features  of  equation  (7a)  are  the  frequency 
indepe dance  of  the  Real  part  and  the  slope  1  dependence  (on  the 
log-log  plot)  in  the  imaginary  part  of  the  impedance  (see  region  I 
on  figure  7). 


(b)  The  intermediate  frequency  range  (wt  »  1  but  r#«T/Rp 


Re(Zo)  =  Im(Zo)  = 


«  1): 
(7b) 
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This  results  in  the  slope  -i  line  in  region  II  of  figure  7. 

(c)  The  high  frequency  range  (mt  »  1  AND  rcwt/R  »  1): 

o  P 

Re(Zo)  =  /Rr 

o 


Im(Zo)  = 


(7c) 


which  results  in  a  frequency  independant  real  impedance,  and  the 
slope  -1  line  in  the  imaginary  impedance  as  can  be  seen  in  region 
III  of  figure  7. 

In  addition  to  the  frequency  limits,  it  was  also  assumed  that 
Vrp  «  1,  and  that  R  >  rg  by  at  least  an  order  of  magnitude  (which 
is  usually  the  case). 

To  combine  the  single  pores  into  a  unified  metal-electrolyte  system, 
we  extend  De  Levie's  model  slightly.  We  define  an  x-y  plane  and 
divide  it  into  Nz  squares,  and  randomly  place  a  rectangular-width  pore 
on  the  x-y  plane  with  its  length  extending  in  the  z  direction.  The 
length  of  the  pore,  L,  is  also  randomly  chosen  to  be  from  1  to  N,  and 
the  pore  width,  W,  is  selected  so  that  there  is  a  Gaussian  distribution 
about  W  =  2  (to  prevent  a  few  very  wide  pores  from  "taking  over”  the 
whole  electrode).  More  random  sized  pores  are  randomly  placed  on  the 
x-y  plane  until  the  desired  porosity  is  reached.  We  may  then  solve  for 
the  total  impedance  of  the  system  by  adding  up  all  the  parallel 
impedances  due  to  the  individual  pores,  using  equation  (5)  to  represent 
the  impedance  of  each  pore.  The  impedances  of  the  non -porous 
sections  of  the  electrode  are  also  added  on  to  the  total  impedance. 

The  results,  for  a  porosity  of  0.3  and  N*40,  are  presented  in  figure 
8.  The  criteria  for  choosing  N  was  based  on  considerations  similar  to 
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those  used  to  choose  N  in  the  Random  Network  model,  namely  that  the 
results  be  independant  of  N  and  statistically  independant  of  the  random 
number.  This  phase  of  the  calculation  was  performed  on  an  IBM  PC. 


III.  DISCUSSION  OF  RESULTS 

There  is  a  qualitative  agreement  between  the  Random  Network  model 
and  the  single-pore  model.  A  comparison  of  the  passive  elements 
representation  of  the  two  systems  (figure  4  and  figure  8)  also  yields 
order-of -magnitude  quantitative  agreement.  One  discrepancy  between 
the  two  models  is  the  definite  shift  In  the  location  of  the  "mid- 
frequency  peak"  in  Figure  8  as  compared  to  the  Random  Network  and 
single-pore  results.  A  major  difference  between  the  two  models  in 
construction,  is  the  fact  that  in  the  Random  Network  the  pores  interact 
with  each  other  as  in  a  true  porous  electrode,  while  in  the  De  Levie 
model  each  pore  is  kept  separate  from  the  others,  and  is  a  separate 
channel.  The  favorable  comparison  of  the  two  models  would  seem  to 
bear  out  De  Levie's  assertion  that  "the  disregard  of  crosslinks  [between 
the  pores]  presumably  does  not  introduce  a  significant  error  [1]. 
Nevertheless,  this  can  be  a  useful  starting  place  for  an  explanation  of  a 
discrepancy  between  the  two. 

The  behavior  of  the  system  impedance  in  the  Random  Network  model 
may  be  somewhat  inferred  from  the  values  of  the  constituent  components 
that  make  up  the  system.  The  low-frequency  (real  and  imaginary) 
Impedance  in  figure  4  reflects  the  values  of  R^  and  C  of  the  unit 
(semiconductor)  interface,  although  the  3008  value  of  Re(Z)  is  noticably 
lower  than  the  IKE  Rp.  In  the  high-frequency  regime,  C  dominates, 
totally  shunting  out  R^,  the  only  contribution  to  the  Real  impedance 
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being  R0 .  This  provides  us  with  an  Indication  that  the  total  behavior 
of  the  porous -metal  system  may  be  inferred  from  the  values  of  the 
components  that  make  up  the  complicated  surface  topography.  Our 
ultimate  goal  is  to  generate  a  detailed  picture  of  the  macroscopic 
interface  from  the  impedance  data  of  the  microscopic  components  that 
make  up  the  system.  This  will  be  considered  further  in  a  subsequent 
paper,  along  with  an  exploration  of  the  physical  significance  of  the 
passive  elements  model  of  the  entire  system. 

We  will  also  examine  the  correlation  between  the  single  pore  analytic 
results  (in  the  various  frequency  limits)  and  the  Random  Network 
model.  Only  a  single  porosity  has  been  considered  thus  far,  but  an 
exploration  of  the  Random  Network  model  at  different  porosities  is  in 
order.  Although  some  rough  surfaces  may  exhibit  a  conduction 
percolation  threshold  [11,12],  we  do  not  expect  to  see  any  classical 
percolation  in  our  system.  This  is  due  to  the  fact  that  we  have 
constrained  the  system  to  be  above  the  percolation  threshold  by  the 
requirement  that  the  metal  electrode  be  continuous.  We  would, 
however,  expect  to  see  some  less  dramatic  system  dependence  on 
porosity. 
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LIST  OF  FIGURES 

Representation  of  the  3-D  Lattice  in  terms  of  Circuit 

elements.  Figure  lb  shows  our  model  in  terms  of  a 

checkerboard  pattern  of  bulk  components  of  metal, 

electrolyte,  semiconductor,  and  air.  Figure  la  depicts  our 

representation  of  (a  corner  of)  Figure  lb  in  terms  of  a 

network  of  circuit  elements.  Z  is  the  impedance  of  the 

sc 

semiconductor  layer  on  the  interface  (see  Figure  2). 
Equivalent  circuits  used  to  represent  the  microscopic 
components  of  the  system  (metal,  electrolyte,  semiconductor, 
etc . ) . 

The  Random  Network  lattice  used  to  generate  the  data  for 
Figure  4.  At  this  porosity  of  .3,  the  electrolyte  does  not 
fully  penetrate  the  lattice.  (Actually,  the  penetration  of  the 
electrolyte  into  the  pores  seems  to  conform  to  classical 
percolation  theory.  At  porosities  above  .35  the  liquid  starts 
to  penetrate  through  the  electrode). 

Impedance  spectrum  from  the  Random  Network  model  at  a 
single  porosity  (0.3),  and  a  single  random  number.  Three 
frequency  regimes  are  identified,  and  we  have  modelled  two 
regimes  with  a  staple  passive  element  model.  The  bars  on 
each  curve  represent  the  variations  in  the  data  when  various 
random  number  seeds  are  used  to  construct  the  lattice.  For 
each  random  number  seed,  the  shape  of  the  Impedance 
spectrum  is  almost  identical. 

The  impedance  in  the  Random  Network  model  as  a  function  of 
system  size.  All  values  are  normalized  to  the  impedance  of 
lattice  size  15  x  15  x  15  (N«15). 


6. 


The  De  Levte  model  of  a  single  pore.  Figure  (6a)  presents 
the  general  transmission-line  model.  (6b)  presents  the  value 
we  have  chosen  for  R,  the  electrolyte  impedance,  and  our 
specific  model  for  Z,  the  electrode -electrolyte  interface  inside 
the  pore. 

7.  Impedance  spectrum  from  the  DeLevie  model  -  single  pore. 
The  sections  of  the  curve  labelled  (1)  are  areas  in  which  the 
Imaginary  impedance  plot  is  of  slope=l  on  the  log-log  plot, 
and  the  impedance  in  these  ranges  may  be  modeled  by  the 
two  passive  elements  models  shown .  The  section  of  the 
curve  labelled  (2)  has  a  slope=l/2;  See  equations  7a  -  7c. 

8.  Impedance  data  from  Dele  vie  model  -  randomly  distributed 
pores.  Porosity  *  0.3,  N=40. 
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